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In this paper we derive, directly from the Nambu-Goto action, the relevant components of the 
acceleration of cosmological featureless p-branes, extending previous analysis based on the field 
theory equations in the thin-brane limit. The component of the acceleration parallel to the velocity 
is at the core of the velocity-dependent one-scale model for the evolution of p-brane networks. We 
use this model to show that, in a decelerating expanding universe in which the p-branes are relevant 
cosmologically, interactions cannot lead to frustration, except for fine-tuned non-relativistic networks 
with a dimensionless curvature parameter fc <g 1. We discuss the implications of our findings for 
the cosmological evolution of p-brane networks. 
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I. INTRODUCTION 

It is generally accepted that the universe underwent 
a phase of accelerated expansion in its early history. 
This paradigm, usually denominated cosmological in- 
flation [l|, [4I, explains the extreme flatness and homo- 
geneity of the observed universe, as well as the the ori- 
gin of the large scale structure. In the context of the 
brane-world realization of string theory, cosmological in- 
flation could be driven by the interaction between p- 
dimensional D-branes 0-13 (which are, along with funda- 
mental strings, the fundamental objects of the theory). 
Branc inflationary scenarios typically end with a symme- 
try breaking phase transition, leading to the production 
of branes with lower dimensionality. The copious pro- 
duction of 1-dimensional branes (cosmic strings) is pre- 
dicted for a large variety of branc inflation models, while 
higher-dimensional p-branes might also be produced [6|- 
|8|]. Therefore, brane inflation may lead to the produc- 
tion of p-branes networks, evolving in higher dimensional 
backgrounds. 

Cosmic strings and domain walls are the only non- 
trivial p-brane solutions allowed in 3 -I- 1 dimensional 
backgrounds. The interest in cosmic strings has recently 
been revived due to the possibility of detecting their cos- 
mological signatures observationally. In particular, they 
may leave an observable imprint on the B-mode polariza- 
tion [^, [lO| and small-scale anisotropy [ll| of the cosmic 
microwave background which may be within reach of the 
Planck mission. Moreover, there is also the prospect of 
the future detection of the gravitational waves emitted by 
cosmic strings with the LIG02 and LISA probes ^12] (see 
also [l3l - [l5| for other examples of possible observational 
signatures of cosmic strings). In order to make precise 
observational predictions it is necessary to understand 
the nature of cosmic strings and their late-time evolu- 
tion. While standard cosmic string networks have been 
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extensively studied both using a semi-analytical velocity- 
dependent one-scale (VOS) model [T^ and numerical 
simulations Il7|, significantly less is known about the dy- 
namics of complex string networks with junctions. For 
example, cosmic superstring networks may have an hi- 
erarchy of tensions and junctions 18* and, consequently, 
their observational signatures are expected to be differ- 
ent from those of ordinary cosmic strings [1^ [20| . Still, 
it is not yet clear whether or not these differences pre- 
vent cosmic superstring networks from attaining a linear 
scaling regime and, therefore, the late-time evolution of 
these networks is not yet completely understood [21h23| . 

A p-brane network, if frozen in comoving coordinates 
(or frustrated), is characterized by a negative average 
pressure. This property is particularly interesting in 
the case of domain wall networks. There is compelling 
evidence [2J, [25| which indicates that the universe is 
currently undergoing a phase of accelerated expansion, 
caused by an exotic energy component (dubbed dark 
energy) accounting for more than two thirds of the en- 
ergy density of the universe and whose nature is yet un- 
known. If domain walls were the dominant energy com- 
ponent, and provided their root-mean-square (RMS) ve- 
locity was small enough, they could drive a phase of ac- 
celerated expansion. For this reason, frustrated domain 
wall networks have been suggested as a dark energy can- 
didate 26] . The conditions under which domain wall net- 
works may frustrate as a result of cosmological evolution 
were intensively studied in j27n-32j using a semi-analytical 
VOS model and field-theory simulations. These studies 
indicate that frustration does not result naturally from 
the evolution of realistic and cosmologically relevant do- 
main wall networks. 

The cosmological dynamics of generic p-brane net- 
works has been studied in [33 using a VOS model char- 
acterizing the evolution of their RMS velocity and char- 
acteristic length. As in the case of cosmic string |3j| 
and domain wall networks [30, |33 , in homogeneous and 
isotropic universes with a decelerating power law expan- 
sion, p-brane networks may evolve towards a linear scal- 
ing solution in which the characteristic length grows pro- 
portionally to the Hubble radius and the RMS velocity 



remains constant [33| . This scale- invariant solution is an 
attractor of the VOS equations, characterizing the late- 
time evolution of the networks in a frictionless regime. 

In this paper we derive the evolution equation for the 
velocity of featureless infinitely-thin p-brane networks in 
homogeneous and isotropic universes with an arbitrary 
number of spatial dimensions, directly from the Nambu- 
Goto action. This generalizes previous results obtained 
using field theory equations in the thin-brane limit, fur- 
ther validating the VOS model for the evolution of p- 
brane networks in iV + 1-dimensional homogeneous and 
isotropic universes derived in [sj. We also consider the 
effect of a generic interaction mechanism between the p- 
branes and other cosmological components, studying its 
potential role in the frustration of the networks. 

Throughout the work, we will assume the metric signa- 
ture [-I-, — , ..., — ] and the calculations will be done using 
units in which c = 1. The Einstein summation conven- 
tion will be used when a greek index variable appears 
twice in a single term, once in an upper (superscript) 
and once in a lower (subscript) position. 



II. FEATURELESS p-BRANES: EQUATION OF 
STATE 



Consider a local inertial frame in which the p-brane is 
instantaneously at rest. If the p-brane is locally flat, a set 
of planar orthogonal coordinates (t,x^ , ■ ■ ■ , cc^) may be 
, x^ parameterize the brane and 



for I ^ i. It is possible to show using similar arguments 
that T°* must also vanish. Moreover, 
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If the p-brane is maximally symmetric with respect to 
the N —p perpendicular directions and its energy is local- 
ized then Derrick's theorem [3a| implies that a necessary 
condition for stability is [33| 
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for ^ > p + 1 (here D = N — p and d^x — dx^'^^ x ... x 
dx^). Spherical symmetry with respect to the D per- 
pendicular directions implies that, at the core, we should 
also have that T'™ — 0, for I >p-l-l, m > p + 1 and 
I ^ m. In most situations of cosmological interest, the 
p-brane thickness is very small compared to its curvature 
radii and may therefore be neglected. Hence, if the p- 
brane is infinitely thin, the non- vanishing components of 
the energy-momentum tensor are 
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where ap is the (constant) p-brane mass per unit p- 

j^ „ . dimensional area, x is a N-vector whose components are 

■■■ X are perpendicular to it. Since the prop- cartesian coordinates, x„ represents the p-brane profile, 



chosen such that 
xP+\ 

erties of a featureless p-brane do not change along the 
parallel directions, there is no component of the physical 
velocity along these directions. Therefore, the compo- 
nents of the energy-momentum tensor T^'^ must be in- 
variant with respect to Lorentz boosts in any direction 
along the brane. 

Consider a boost along one of the parallel directions 
x^ {i — I,--- ,p throughout the paper). The energy- 
momentum tensor, transforms as 



and (5^(x) is the TV-dimensional Dirac delta function. 

The equation-of-state parameter of a brane gas can be 
obtained by performing a Lorentz boost to the energy- 
momentum tensor of a single static flat p-brane and, 
then, averaging over all possible orientations of the brane 
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where 
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A;!-A^=7, A?=A^-7«, 



A; = 1 , for I ^ i , 

and all other components vanish. Here / = 1, ■ 
Hence 

T°''' = 7r°' + jvT'-^ = r"' , 

T*''' = 7T'' -I- 7wT°' = r*' , 
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which leads to 
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Here, v is the RMS velocity of the branes, p = V~^ J pdV 
is the average brane density, P — V~^ J T^dV is the av- 
erage brane pressure and T^ is a large volume. 

Note that Eq. ([T3)) has two important limits. In the 
relativistic limit, with u — > 1, one has w — > 1/iV inde- 
pendently of the value of p. In the non-relativistic limit 
{v — > 0) one has w = —p/N. Taking into account that 
the Raychaudhury in a A''-|-l dimensional FRW universe 
is given by 
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where pi, is the background density, Vb is the background 
pressure and Gn+i is the A^ + 1 dimensional Newton 
constant, we conclude that only domain walls (p-branes 
with p = A^ — 1) or a cosmological constant (p-brane 
with p = N) could lead to an accelerated expansion of 
the universe. 



III. p-BRANE DYNAMICS 

In the zero-thickness hmit, the world- history in space- 
time of a featureless p-brane may be represented by 



x^' =x^'{u''), 



(15) 



where u'^ with D — 0,1, ...,p are the coordinates parame- 
terizing the (p -\- l)-dimensional worldsheet swept by the 
p-brane, u'^ is a timelike parameter and u* are spacelike 
parameters. The action of the p-brane may be derived 
from the underlying field theory 
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where g = det(g^y), g^i, is the metric tensor and C is the 
lagrangian of the model. In the vicinity of the p-brane, 
the line element can be written as [391 



ds^ — gjj_iidx'-'dx'^ + dr ■ dr , 
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is the worldsheet metric with 



x"^ — dx"" /dufi and (dr • dr)^/^ is the infinitesimal dis- 
tance to the brane (p-dimensional) core. Therefore, the 
volume element is given by 
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Since one is assuming that the p-brane is thin and fea- 
tureless, the lagrangian density can only vary along the 
perpendicular directions and, as consequence, it depends 



only on the cc^"*"^, • • ■ 



,-N 



coordinates. Integrating the 



action with respect to these coordinates, one obtains the 
Nambu-Goto action for infinitely thin p-branes 
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where 
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is the (constant) p-brane mass per unit p-dimensional 
area. 



IV. EQUATION OF MOTION 

By varying the Nambu-Goto action in Eq. (fTOj) with 
respect to x^^, one obtains the equations of motion for 
the dynamic variables a;^ 
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In a {N + l)-dimensional fiat Friedmann- Robertson- 
Walker Universe, the line element is given by 



ds^ = a^{rj) {d^rj — dx ■ dx) , 



(22) 



where a is the scale factor, i] — J dt/a is the conformal 
time, t is the physical time and x is a N- vector whose 
components are comoving cartesian coordinates. Since 
the Nambu-Goto action is invariant under worldsheet pa- 
rameterizations, one is free to impose temporal transverse 
gauge conditions 
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A dot represents a derivative with respect to conformal 
time and x^ = 9x/9mj. These gauge conditions are cho- 
sen so that the timelike worldsheet coordinate is identi- 
fied with the conformal time and x represents the phys- 
ical velocity, perpendicular to the p-brane itself. More- 
over, the local orthogonal coordinate system is chosen in 
such a way that the coordinate lines coincide with the 
principal directions of curvature. Eq. (|2ip yields 
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and % = a/ a. 

Let us define the unitary vectors 
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where w = |x| is the velocity of the brane at a particular 
point and e^ (with i = 1, ...,p) are the unitary tangent 
vectors along the principal directions of curvature. Dif- 
ferentiating Eq. (j24l) with respect to conformal time one 
finds 

a^ = (x • e.) ej = w je^ , (29) 



and, consequently, there is a component of the accelera- 
tion parallel to the p-brane 



Using Eq. (|25l) one may show that 
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where we introduced the comoving curvature vector along 
the principal direction of curvature i defined by 
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(32) 



n::— constant 



with j ^ i- Here we introduced the physical length along 
the i direction ds^ = x^ du^. Using Eq. pi I) , taking into 

account that x = wv + wv and that v is perpendicular to 
V, one finds that the tangential acceleration (parallel to 
the velocity) is given by 
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where we defined the tangential curvature as the projec- 
tion of the comoving curvature vectors along the velocity 
direction kj« —kq-w. This tangential acceleration allows 
us to obtain an evolution equation for the velocity of the 
p-brane 



+ {l-v^) [(p + l)Hw-K||] =0, 



(34) 



where the total tangential curvature, kii = X^^-i ^iih ^^^ 
introduced. This equation is identical to that obtained 
from field theory equations in Ref. [33 • 

There are N — p — 1 directions which are simultane- 
ously perpendicular to the p-brane and to its velocity. 
Let us denote the unitary vectors along these directions 
by e; with I = p + 2, ..., N. The acceleration along these 
directions is given by 

p 
a; = {x ■ ei) ei = ei{l - v'^)^k-i ■ ei = 

= ei{l~v^)K^i, (35) 

where we have introduced the total comoving curvature 
along the perpendicular direction I k±i — ^~_,ki ■ ej. 
Therefore, the total perpendicular acceleration is given 
by 
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with X = x= -|- X|| -|- xx 



V. VOS MODEL 

Let us consider a network of p-branes in a {N -1-1)- 
dimensional FRW Universe and define the root mean 
square (RMS) velocity, v = \/{v'^) as 
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where d^u ~ du^ x ... x du^ . The characteristic length, 
L, of the network is defined as 



(38) 



where p is the average brane density. 

A unified VOS model for the dynamics of p- brane net- 
works in (A^-l-l)-dimensional FRW universes was derived 
in [33|. This model is described by the following equa- 
tions 



dv 
'di 
dL 

m 



(1- 
HL 



V 



L 

DL 



-V +-V, 



(39) 
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where £'^^ = {p + l)H is the damping lengthscale, H — 
Ha is the Hubble parameter, D = N — p, c > is the 
energy-loss parameter, k — R\\L/a is a dimensionless cur- 
vature parameter. 
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and the assumption that («■*) — w^ was made (see 
|34|V A frictional force — caused by the interaction of 
the branes with ultrarelativistic particles or other fric- 
tional sources — may be included in Eq. (|5n| . by in- 
troducing an extra term in the damping lenghtscale, 

£^^ = {p + i)H + ej\ 

Consider the presence of an interaction mechanism be- 
tween the p-branes and a component average density pmt ■ 
A very conservative upper limit to the total momentum 
per unit volume transfered from that component to the 
p-branes in one Hubble time. 
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is given by pint- In this case one has 

<XH, 



(42) 
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where x = Pint/p- 

li a (X t^ (with < /? < 1), and the friction scale is 
negligible compared to the Hubble radius (£/ ^ H~^), 
this model may admit linear scaling solutions of the form 
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See [33 
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for a discussion of other scaling regimes. 
= constant, then 
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The effects of this interaction mechanism may slow down 
the branes slightly but, if the interacting component of 
energy density pint is subdominant (that is x < 1), its 
potential role on the frustration of the network (defined 



by L <C H~^ and u ^ 1) is very limited. Relaxing the 
assumption that the p-brane network is the dominant 
energy component may help frustration: if x 3> 1 (or 
equivalently if pint 3> p) the interaction mechanism may 
decelerate the branes effectively, leading to the frustra- 
tion of the network. However, in this case, pint would be 
the main contributor to the energy budget. Of course, 
assuming that x is time independent is unrealistic, since 
one would expect the expansion of the background to af- 
fect the efficiency of any realistic interaction mechanism 
(in particular, the ratio x = Pint/p is expected to be a 
function of a). Nonetheless, Eq. ([M)) shows that consid- 
ering a time varying x does not help much if its present 
value Xo ^ 1- As a matter of fact, frustration might only 
occur, under these circumstances, for networks which 
have fc <C 1 for w <g; 1. This gene ralizes a well known 
result for domain walls [23, [3l|, [S^l to arbitrary p-brane 
networks. In the particular case of domain walls [3l|, 
there is very strong analytical and numerical evidence 
that domain wall networks (with or without junctions) 
arc unlikely to attain fc <IC 1 for w ^ 1, if they are the 
dominant energy component. This effectively rules out 
domain walls as a cosmologically relevant dark energy 
candidate: frustration can only occur either if the net- 
work is designed to have fc ^ 1 in the non-relativistic 
limit — which appears to be unrealistic — or if x is 
much larger than unity — in which case the domain wall 
energy density would be subdominant. For p < N — 1, 
the damping is less efficient, and thus frustration is even 
less likely to result from the natural evolution of the net- 
work. Therefore, unless there was a natural mechanism 
that drives k towards zero in the non-relativistic limit 
(which seems unlikely) the "no-frustration conjecture" is 
also expected to apply to any realistic and cosmologically 
relevant p-brane network. 

In ref. [40| the authors perform field theory simulations 
of a model with Z2 x U{1) symmetry in (2+l)-dimcnsions. 
Their model has two discrete vacua, allowing for domain 
walls and a conserved Noether charge. The authors argue 
that the Noether charge and currents become localized on 
the walls, forming kinky vortons and providing a possible 
mechanism for the frustration of domain wall networks. 



However, the authors never calculate the overall equation 
of state of the network. Had they done that, they would 
have found significant deviations with respect to that of 
a frustrated featureless domain wall gas (w = —2/3). 



VI. CONCLUSIONS 

In this paper, we derived the equation of motion for in- 
finitely thin featureless p-branes, by computing the tan- 
gential and normal components of the acceleration di- 
rectly from the Nambu-Goto action. Our results further 
validate the semi-analytical VOS model developed in [33 , 
and its use in dynamical studies of p-brane networks. The 
VOS model unifies, in a single framework, the dynamics 
of p-brane networks for any possible values of the pair 
{N,p). While part of the dynamical dependency on the 
parameters N and p is explicit in the equations of motion, 
there is also an implicit dependency in the parameters c, 
£f and k. In this paper we demonstrated that, if the 
p-branes are the dominant component of the universe, 
then frustration is not possible except if the curvature 
parameter is driven towards very small values for non- 
relativistic networks or if the expansion is accelerated. 
In the case of domain walls there is very strong analyti- 
cal and numerical evidence (both in two {N — 2,p = 1) 
and three {N = 3,p = 2) spatial dimensions) that k is 
never becomes much smaller than unity (except deep into 
inflationary or friction dominated regimes), thus prevent- 
ing frustration from being attained, at least for realistic 
domain wall networks playing a dark energy role. We 
conjecture that this may be a general result, valid for 
any realistic p-brane network independently of the val- 
ues of N and p with 1 < p < N — 1. 
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